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ABSTRACT 

An explicit example is presented (a one-loop triangle graph) where dimensional regular- 
ization fails to regulate the infra-red singularities that emerge at intermediate steps of 
studying large-Q 2 Sudakov factorization. The mathematical nature of the phenomenon 
is explained within the framework of the theory of the As-operation. 
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1 Introduction. The power of dimensional regularization in applied Quantum 
Field Theory is nothing short of miraculous. Its notational and ideological economy makes 
it a perfect tool for a wide range of perturbative calculations^ as well as a convenient 
instrument for theoretical studies of factorization problems^. 

Practical calculations in perturbative QCD typically exploit the presence of a large 
kinematic variable Q and deal with (the leading terms of) the corresponding asymptotic 
expansion of the cross section. Results describing the leading terms of such expansions 
generalize the familiar Wilson short-distance operator product expansion and are known 
as factorization theorems (for a review see e.g. ||). 

The starting point of the construction of such an asymptotic expansion is a formal Tay- 
lor expansion of the integrand with respect to the asymptotic parameter of the problem. 
Such an expansion typically generates infrared singularities, in addition to the ultraviolet 
divergences of the original diagrams. Even though there exist well-defined prescriptions 
for eliminating the infrared singularities in the (Wilson) coefficient functions of an expan- 
sion, the singularities are important at intermediate stages. 

The main advantage of dimensional regularization is that it simultaneously regulates 
both ultraviolet and infrared divergences while preserving Lorentz and gauge invariance. 
Situations where dimensional regularization works well comprise all expansion problems of 
a Euclidean type [f|,0 and many problems of an inherently non-Euclidean (Minkowskian) 
nature 0. 

The aim of the present note is to show that there are singularities that arise in the 
asymptotic analysis of Minkowski space problems and that cannot be regulated by dimen- 
sional continuation. We will pinpoint the origin of the difficulties by means of an explicit 
one-loop example: a form factor graph. We will show that the nature of the problem 
is very general so that the same difficulty is bound to arise in many Minkowski space 
situations. 

A subsidiary aim is to popularize the concepts of a new and very general method of 
analyzing asymptotic behavior of Feynman graphs — the method of the /ls-operation for 
products of singular functions |jj . The theory of the As-operation has been worked out in 
detail for the Euclidean problems J7|, ||, where it provides a simple and clean method 
for treating such problems as renormalization, the operator product expansion, and the 
large-mass expansion. While the principles of As-operation |7j are very general, there are 
several technical problems that one has to overcome before a fully satisfactory extension 
of the method to the most general non-Euclidean situations is achieved. The results we 
describe represent a first step in that direction. 



2 The example. Consider the one- loop diagram of Fig. 1. We will see that neither 
the possible presence of UV divergences nor the structure of the numerator of the cor- 

1 See, e.g., the record-setting 5-loop calculations in the <p% =4 model the multiloop calculations in 
QCD [Q; the industrial-scale QCD calculations of multijet hadronic scattering processes ||. 

2 See, e.g., the derivation of the general Euclidean asymptotic expansions of perturbative Green func- 
tions §,!,§. 
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responding integrand are relevant; only the structure of the denominators is important. 
Therefore, it is irrelevant whether the corresponding particles are scalars, spinors or vec- 
tors. For concreteness, one can imagine that the horizontal fat line is a gluon in Feynman 
gauge while the side lines are quarks, and we will use these names for the lines. We 
consider the asymptotic behavior when Q 2 = — q 2 = — (pi — P2) 2 gets large with all other 
invariants — p 2 , p\ and the masses — held fixed. 

To focus on the analytical effect we wish to describe, we will choose the quarks to be 
massless, and the gluon to have mass m,Q although the choice of masses does not affect 
the general principles of our method. The external quarks are on-shell: p\ = p\ = 0. 



We assume that the MS scheme |TT| is used for UV renormalization; the dependence 
on the renormalization parameter \x thus introduced is logarithmic and known explicitly 
(the renormalized diagram has the form of //-independent terms plus constant x log/z). 
Anyhow, since we will eventually concentrate on studying the integrand prior to loop 
integration, the UV behavior is also inessential. 

There are two essential dimensional parameters in the problem: Q 2 and m 2 . The third 
parameter is the renormalization scale /1. The Sudakov asymptotic regime is Q 2 — > 00, 
with fixed m 2 and /i. By dimensional analysis this is equivalent to m 2 ,/x 2 — > with fixed 
Q 2 . (For definiteness, we will consider the case of a space-like momentum transfer q^.) 
Since the dependence on /1 is known explicitly, it is sufficient to consider the expansion 
at m 2 — > with Q 2 and /1 fixed (cf. [|13[). This is the most convenient way to proceed 
within the techniques of the As-operation |7j. 

The integrand for the graph is 

111 

l(yki''P\ P2 777 j —— X X (-0 

k 2 — m? + ii] {k — pi) 2 + it] (k — P2) 2 + ii] 

Considering contributions to the integral from various regions of integration space is 
equivalent to considering various integrals of the form 

G[(/>;p 1 ,p2,m] = J d D kI(k;p 1 ,p 2 ,m)(f)(k), (2) 

where <fi is an arbitrary test function (i.e. a smooth function which is non-zero only in 
a finite subregion of the integration space). So, one arrives at the basic problem of 
expanding arbitrary integrals of the form (0) in m — > with pi and p 2 fixed and light- 
like. This is exactly the same as to say that one has to expand the integrand ([!]) in 
powers and logarithms of the small parameter m 2 in the sense of distributions. Since 
such expansions commute with multiplication by polynomials 0, one can forget about 
numerators of propagators as well as possible vertex factors. 

Construction of an expansion in the sense of distributions involves the following steps 
|7j]: formal Taylor expansion of the product (|TJ) in powers of m; classification of sin- 
gularities of such a formal expansion; construction of counterterms to be added to the 

3 Note that a non-zero gluon mass may have a dynamical origin |fl0| . Gauge invariance, however 
important for studying combinatorial properties of entire perturbation series, is of no relevance to the 
analytical problem of asymptotic expansions of individual graphs. If one wishes, one can consider the 
graph to be one in an abelian gauge theory, where the gluon can consistently be given a mass. 
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formal expansion to transform it into a correct expansion in the sense of distributions. 
Construction of counterterms proceeds in an iterative fashion: from singularities of sim- 
pler nature — i.e., to which fewer singular factors contribute — to the more complicated 
ones; construction of counterterms for the latter involves simpler counterterms obtained 
at previous iterations. 

For practical reasons it is convenient, whenever possible, to employ a regularization 
throughout the entire procedure. Dimensional regularization is the prime candidate for 
that role. However, as we will see shortly, it fails to regulate a singular expression at an 
intermediate step of constructing the expansion in the sense of distributions for (p. 

3 Geometry of singularities of the formal expansion. We follow the general 
procedure that gives the As-operation. First we make a formal expansion of the integrand 
in powers of m. This gives a correct leading order expansion of the integral with a test 
function, provided that the test function is zero in a neighborhood of all singularities of 
the expanded integrand. That is 

f 1 1 1 

G[(p;p 1 ,p 2 ,m} = / d D k(p(k) x — — x , 

J k -\- ITj [k — Pi) +M] {k — P2) + IT) 

+ o(l), if 4> is zero near singularities. (3) 

Here the remainder is a power of m smaller than the leading term. 

To get a complete expansion, valid for all test functions, one should examine a small 
neighborhood of the singularities of the first term on the right of Eq.(|]). The principles 
of the As-operation prescribe that we should start with the highest dimension singular 
surface and then treat successively lower dimension surfaces. 

Since the integration contours may be deformed away from the light-cone singularity 
of any single propagator (except at its apex), such a singularity is integrable. Hence the 
only singularities that we need be concerned with are the apexes of the light cones and 
the intersections of the light-cone singularities of different propagators, 

Now, the first denominator in Eq.(^]) generates a singularity localized on the light-cone 
k 2 = 0. The second factor is singular on (k — pi) 2 = which is nothing but the light-cone 
k 2 = shifted so that its apex is at the point k = p\. Since p\ is a light-like vector, the 
two light cones intersect on the straight line 

A = {k = zxpi, —oo < z\ < +oo}. (4) 

This intersection is non-trivial because the light cones are not transverse at the intersection 
points. Similarly, the singularities of the first and third factor overlap on the line 

B = {k = z 2 p 2 , -00 < z 2 < +00}. (5) 

The singularities of the second and third factors intersect on a smooth manifold and 
are transverse there; this implies that in a small neighborhood of each intersection point 
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the integral factorizes, so that the singularity is integrable and such an intersection is 
harmless. The exception is the point k = where the first denominator is also zero, and 
we consider this point separately. 

There are three points — namely, k — 0, k — pi, and k = p2 — where the analytical 
nature of singularities is more complicated than at the generic points of the lines 
In particular, the singularity 

S = {k = 0}, (6) 
is where the effect we are after takes place. 

The geometrical pattern of singularities can be visualized as in Fig. 2. 



4 Structure of singularities at k oc p\. Consider singularities near a generic point 
on the line A, (^). Fix z ^ 0, 1, and consider a small neighborhood O of the point zp\. The 
third factor of the integrand is smooth in O and, therefore, can be effectively relegated 
to the test function. So, within O it is sufficient to study the expansion of the product 
of only the first two factors that contribute non-trivially to the singularity in O. It is 
convenient to choose light-cone coordinates, k = (k + ,k_,k±), so that p\ = (pi+, 0, 0j_), 
P2 = (0,p2~, 0j_). Our conventions will be such that k 2 = k + k- — kj_, 2k -p\ = pi + k_, and 
d D k = \dk + dk_d D ~ 2 k±. 

Then the formal expansion of the first two factors takes the form: 

1 1 



z Pi+k~ — k\ — m? + ir] (z — l)pi + k_ — k\ + ir] 

1 1 

= ) T2 — - x t ^ ; T2 — - + ( 7 ) 

zp\ + k^ — k^+irj (z — ±)pi + k- — k L + ir] 

One can see that the singularities of the product on the r.h.s. are localized at the origin 
of the space of the variables k_ and k±. This means that ([?]) holds in the sense of 
distributions on test functions that are zero near = k± = 0. The situation here is very 
similar to what one has in the case of a single Euclidean propagator treated in section 
7 of @ — see especially eqs. (7.25) and (7.26). The only important difference is that the 
singular functions on the r.h.s. are not homogeneous. This is remedied by the change of 
variable k_ = ±t 2 after which the expressions on the r.h.s. of (0) become homogeneous 
with respect to simultaneous scaling in t and k±, and simple power counting shows that 
the singularity is logarithmic (at D = 4). 

This allows one to repeat the reasoning of section 7 of 0] to write down the following 
analogue of eq. (7.25) of 0: 

The l.h.s. of eq.(7) = 



zpi + k^ — k 2 i + ir] (z — l)pi + k_ — k 2 i + ir] 
+ c A (m 2 , z )—S(k^ D - 2 \k ± ) + o(l), (8) 



Pi- 
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with 



CA(m 2 ,z) = pi + J dk^d D 2 k A 



zpi + k_ — k\ — m 2 + it] (z — l)pi + A;_ — k\ + if] 
1 1 



Zpi + k^ — k 2 L + il] (z — l)j9i + /c — k\ + 17] 
r 1 1 

J ± zpi + k_ — k\ — m? +17} (z — l)pi + k- — k\ + irj ^ 

In the second form of this integral we have used the fact that the integral of a homogeneous 
function^ is zero within dimensional regularization. 

At this point dimensional regularization works. Recall [0] that the role of the coun- 
terterm ca is two-fold: first, it cancels the divergence of the first singular expression (after 
integration with test functions); second, it ensures that the expansion is asymptotic (o(l)) 
on any test function. The expansion holds in the sense of distributions on any test 

function (in the variables k±). The simplification gained is that, whereas the l.h.s. of 
(§) has a complicated dependence on m, the first two terms on the right have a simple 
power dependence, which becomes logarithmic at D = 4. This we will verify from the 
explicit calculation of Ca in the next section. 



5 Explicit expressions for counterterms. It is not difficult to perform the inte- 
grations in (|D explicitly^ to obtain 

c A (m 2 , z) = 2i6(<d<z< 1) T(e) tt 2 ^ [(1 - z)m 2 ] , (10) 

where e = (4 - D) /2. 

There are a few points worth making here. First, the counterterm is zero for z > 1 and 
z < 0. This agrees with the fact that according to the Landau equations^] the singularities 
at those values of z are not pinched. In the language of distribution theory, this says that 
for z < and z > 1 the product in the first term of the r.h.s. (Bl) is (i) well-defined in the 



4 Recall that the reasoning that leads to this expression is, strictly speaking, done in the "straightened" 
coordinates (t,k±) with homogeneous functions. 

5 Recall that the definition of dimensional regularization ffl prescribes that we should perform the 
integrations over the transverse components first. In the present case, however, the same result is re- 
produced in a slightly easier way if one first performs integration over k- . At this point one should be 
aware of the fact that the underlying definition of the integral is in terms of the homogeneous coordi- 
nates, and the cutoffs that are present at intermediate stages (cf. @]) are symmetric with respect to the 
homogeneous coordinates and become asymmetric in the coordinates fc_ , k± . Such subtleties stress the 
need for a meaningful rigorous definition of dimensional regularization in momentum/coordinate space 
representations. Some of the issues involved will be discussed in a forthcoming publication 

6 See e.g. ^| . Note that the Landau equations are usually associated with studying analyticity 
properties of Feynman diagrams. Finding cuts etc., however, is equivalent to determining when an 
expansion near the corresponding value of external momenta and masses contains non-trivial (i.e. non- 
analytic) contributions. In the context of the expansion problem proper this issue was reconsidered by 
Libby and Sterman [ fl3[ . A reinterpretation of the Landau equations from the point of view of asymptotic 
expansions of distributions is presented in E3 . 
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sense of distributions and (ii) does not require additional counterterms to approximate 
the l.h.s. in the sense of distributions to o(l). (In general, (ii) does not necessarily follow 
from (i).) 

Second, the counterterm (|T0|) is the only term in (H) that contains a non-analytic de- 
pendence on m (cf. the discussion of the role of such counterterms in 0]). After expansion 
in e the dependence becomes logarithmic .[] 

Third, the expression ( |10D and all its derivatives in z have well-defined limits as 
z — > +0. This conclusion does not change if the second (quark) propagator in (|I|) and, 
correspondingly, in the first term of (^) contains a non-zero mass, say, mi, because then 
the last factor in flT0| ) will simply be replaced by [zm 2 + (1 — z)m 2 ] t . Note, however, that 
the presence of mass in the first factor (the bottom line or gluon in Fig. 1) is important. 

Finally, it is not difficult to rewrite the above expansion in an explicitly covariant 
form: 

1 1 

V 

k 2 — m 2 +irj (k — pi) 2 + irj 

= Vrr-^7Z 1IX -+ C dz5 (D \k-z Pl )c A (m 2 ,z)+o(l). (11) 

k + ITj (k — pi) z + irj Jo 

We have explicitly taken into account that ca vanishes outside the interval < z < 1. 
Note that the 5-function on the r.h.s. is D-dimensional. This expression represents the 
leading power term in a correct asymptotic expansion in the sense of distributions on test 
functions that are zero in small neighborhoods of k = and k — pi. Note that a similar 
expansion is obtained for the expansion of the product of the first and third factors of ([!]) 
(with pi replaced by p 2 ). 



6 Taking into account the third factor. Let us now consider the entire expression 
([I]). In order to transform its formal expansion, the integrand of (|3|), into a well-defined 
expansion in the sense of distributions, the general recipe of the Extension Principle of 
tells to add counterterms localized at singular points of the formal expansion in (||]) 
with properly chosen coefficients. Let us show that the expansion that is valid on the test 
functions that vanish in neighborhoods of the points k — 0, k — pi and k = p 2 is given 
by the following formula: 

1 1 1 

v x 

k 2 — m 2 + irj (k — pi) 2 + if] (k — p 2 ) 2 + if] 
1 1 1 

— v x 

k 2 + it] (k — pi) 2 + it] (k — p 2 ) 2 + if] 

1 f 1 

+ n \2 I • x / dz 5 {D \k - zpx) c A (m 2 , z) 

{k — p 2 ) z + ti] Jo 

+ 77 \tT~ x f 1 dz5( D \k-zp 2 )c B (m 2 ,z) + o(l). (12) 

{k — piY + iri Jo 

7 It is possible to avoid the use of dimensional regularization in expressions like — cf. |J where 

the results of H are presented in a regularization-independent form. 
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Indeed, any test function <p(k) that vanishes in small neighborhoods of the points k = 0, 
k = pi and k = p 2 , can be represented as a sum (fx + <p 2 where <pi is zero around the 
line k oc P2 and tp 2 is zero around the line k oc p%. On (fx, the second counterterm 
vanishes and we are left with the product of the expansion ( |TT1) times the third factor. 
Since the product of the third factor and <pi is a valid test function <px, one arrives at a 
correct expansion to o(l). A similar reasoning is applied to (p 2 (note that the function cb 
coincides with ca in our example). 

It follows that the expansion ([12]) is actually correct for all test functions that vanish 
in small neighborhoods of the points k = 0, k = p\ and k = p 2 . To make the expansion 
valid on all test functions, one has to add appropriate counterterms localized at the points 
k = 0, k = pi, and k = p 2 . Such counterterms are, in general, linear combinations of 
(^-functions and their derivatives localized at those points with coefficients depending on 
the expansion parameter in a non-analytic (logarithmic) manner [^]. 

To construct the additional counterterms, one must (i) perform an appropriate power 
counting in order to determine the strength of the singularity; (ii) introduce an interme- 
diate regularization to make the singularity manageable (in Euclidean problems J7| di- 
mensional regularization automatically regulates all singularities) or perform an explicit 
subtraction (as in 0); (Hi) determine an explicit form of the counterterms that need to 
be added in order to ensure the approximation property of the resulting expansion. 

It will be at step (ii) that dimensional regularization fails in our example. 



7 Singularity at k = 0. Let us focus on the point k = 0. One has to study the 
singularity of the entire r.h.s. of (0) at k — > 0. The r.h.s. of (|1^) contains contributions 
that are analytic in m, and those that are not. It is sufficient to consider the latter 
since they cannot be affected by how the analytic contributions are treated, and it is the 
non-analytic terms that will exhibit the failure of dimensional regularization. 

The terms with non-analytic dependence in m are known explicitly in our case: 



1 



x / dz 5 {D) (k - zp x ) c A (m 2 , z) + (1 <-> 2) 



(k - p 2 ) 2 + it] Jo 

i C dz 5 {D) (k - zp x ) c A (m 2 , *) + (!«-> 2) 



Q 2 Jo z — ir\ 
One can immediately see that: 



(13) 



(a) One has to deal with the product of a one-dimensional distribution l/(z — ii]), 
which itself is well-defined if integrated with smooth test functions, times 0(z). The 
resulting expression is singular and ill-defined at z — 0. The distribution l/(z — ir]) 
is generated from the propagator l/[(k — p 2 ) 2 + irj] when we set k = zp\. 

(b) Dimensional regularization does not regulate this singularity because the form of 
the product is independent of e and there remain no "unused" extra dimensions. 
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(c) There are no cancellations at k = between contributions from the two counterterms 
corresponding to the two singular lines. 

It remains to note that the effect of breakdown of dimensional regularization in the 
above example persists (even if the formulas become more cumbersome) if one introduces 
masses of order m — say, mi and m>2 — into the two quark propagators, or allows the 
external quarks to be off-shell by 0(m 2 ). One can also see that the same configuration 
of singularities emerges e.g. in the studies of the large- s limit (see Fig. 3). All this points 
to universality of the phenomenon of the breakdown of dimensional regularization in 
expansion problems in non-Euclidean asymptotic regimes. 

8 Conclusions. We have considered a rather typical non-Euclidean expansion prob- 
lem (a one-loop form factor graph in the Sudakov asymptotic regime) within the frame- 
work of the theory of As-operation, and we saw some significant differences from Euclidean 
problems. In particular, we have discovered a class of singularities which are not regulated 
by dimensional regularization. 

The origin of the dimensionally-nonregularizable singularities in non- Euclidean asymp- 
totic expansion problems is completely general: 

First, the non-trivial ("pinched") singularities of the expanded integrands in the case 
of non-Euclidean asymptotic regimes may be localized on manifolds with boundaries 
(which is never the case for Euclidean regimes, where singularities are always localized on 
linear subspaces of the space of integration momenta). 

Second, construction of a complete expansion requires introduction of counterterms 
that contain the non-analytic dependence on the expansion parameter and are localized 
on such manifolds; such counterterms may have coefficients that, together with their 
derivatives, possess finite non-zero limiting values at the boundaries of such manifolds 
when the boundaries are approached from within the manifold (z — > +0 in our case) 
while being identically zero outside the boundary. In other words, if the boundary is 
described by the equation z = in local coordinates with z > corresponding to the 
pinched submanifold, then the coefficients near the boundary have the form 8(z)x (a 
smooth function of z) — even for D^4. 

Third, light-cone singularities of the factors that do not contribute to such countert- 
erms may pass over the boundaries of the corresponding manifolds; when projected onto 
such manifolds (which is exactly what happens when one introduces the counterterms 
into the entire product — cf. (|l"3"D), they take the form l/(z ± irj). 

Fourth, such manifolds may be geometrically positioned so that the extra dimensions 
that are instrumental in the mechanism of dimensional regularization are "used up" in 
the counterterms and do not provide any suppression for the resulting singular product 
of the type 8(z) x l/(z±irj). 

It should be emphasized that the problem here is not an ambiguity as in the case 
of 75, but a failure of dimensional regularization to regulate a particular class of infra- 
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red singularities. Moreover, since expansions in the sense of distributions in powers and 
logarithms of the expansion parameter are unique, one has to conclude that it is impossible 
to get rid of the problem by choosing a different "factorization scheme" . 

The problem is certainly associated with our insistence on strict power-and-logarithm 
expansions: The 1/z singularity gives a problem because we have expanded everything 
else in the integrand in powers of a small variable. However, the requirement that the ex- 
pansions to be constructed run in powers and logarithms of the expansion parameter (the 
requirement of "perfect factorization" 0) cannot be relaxed for both phenomenological 
and technical reasons. In particular, such expansions possess the property of uniqueness 
which greatly facilitates iterative construction of the expansions, relieving one of having 
to worry about unitarity and gauge invariance of the final results etc. [[/J. This is partic- 
ularly true in non-Euclidean problems, since it is only at the leading twist level that we 
get simple factorization theorems.^ 

On the other hand, the existence and nature of the anomalies one may have to deal 
with as a result of the effect we have described is not obvious. One thing is clear, how- 
ever: whether one opts for other regularizations (e.g. analytic regularization which would 
replace l/(z±ii]) in the above expressions by (z ± ii])~ x ), or chooses to combine di- 
mensional regularization with a formalism involving direct subtractions as in || — or to 
forgo dimensional regularization altogether in favor of the latter — the consequences may 
be rather unpleasant both for practical calculations and for the general theory of higher- 
twist factorization. 
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8 In the method one of us has given [|16| for treating the Sudakov form factor at the complete leading 
twist level, non-dimensionally-regularizable singularities are avoided either by the use of axial gauge or 
by the use of an equivalent trick in Feynman gauge. One consequence of the resulting lack of "perfect 
factorization" is an annoying proliferation of remainder terms in Ward identities. These are especially 
tricky to handle in a non-abelian theory. 
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Figure captions. 

Fig. 1. The triangle graph corresponding to our example ([!]). The bottom line corresponds 
to a gluon with non-zero mass m, and the side lines to massless quarks. 

Fig. 2. The geometrical pattern of singularities due to the three denominators of the formal 
expansion, (H). 

Fig. 3. A configuration of propagators (and singularities) essentially similar to that in Fig. 1 
emerges in the large-s small-t problem. 



